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Abstract 

We use an effective field theory (EFT) which contains only short-range interactions to study the 
dependence of a variety of three-nucleon observables on the pion mass. The pion-mass dependence 
of input quantities in our "pionless" EFT is obtained from a recent chiral EFT calculation. To the 
order we work at, these quantities are the 1 So scattering length and effective range, the deuteron 
binding energy, the 3 Si effective range, and the binding energy of one three-nucleon bound state. 
The chiral EFT input we use has the inverse 3 Si and 1 So scattering lengths vanishing at m c ^ %t = 
197.8577 MeV. At this "critical" pion mass, the triton has infinitely many excited states with an 
accumulation point at the three-nucleon threshold. We compute the binding energies of these states 
up to next-to-next-to-leading order in the pionless EFT and study the convergence pattern of the 
EFT in the vicinity of the critical pion mass. Furthermore, we use the pionless EFT to predict 
how doublet and quartet nd scattering lengths depend on m n in the region between the physical 
pion mass and m n = m'^" lt . 
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I. INTRODUCTION 



The NPLQCD collaboration recently computed NN correlation functions in QCD by 
using Monte Carlo techniques to evaluate the QCD path integral on a discrete Euclidean 
space-time lattice [l|. This provided the first calculation of nuclear physics quantities from 
full QCD: the NN scattering lengths in the 3 Si and 1 So channel. However, the NPLQCD 
computation was performed at quark masses which are significantly larger than those in the 
physical world. This leaves us with the challenge of understanding how observables such as 
the NN scattering lengths depend on parameters of QCD such as m u and rrid- 

This challenge can be addressed using effective field theories (EFTs). EFTs allow the 
calculation of physical observables as an expansion in a small parameter that is a ratio 
of physical scales. In this paper, we will discuss two different EFTs. The first describes 
the low-energy sector of QCD by exploiting QCD's approximate chiral symmetry This 
chiral EFT (%EFT) is a powerful tool for analyzing the properties of hadronic systems at 
low energies in a systematic and model-independent way. It is formulated in an expansion 
around the chiral limit of QCD which governs low-energy hadron structure and dynamics. 
In xEFTs, the quark-mass dependence of operators in the effective Lagrangian is included 
explicitly. Loops also generate non-analytic dependence on the quark mass, but their effects 
can be computed in a controlled way, since they can be obtained up to a given order in the 
expansion parameter m q /A x , where A x ~ m p is the scale of chiral symmetry breaking in 
QCD. Over the past 15 years, considerable progress has been made in understanding the 
structure of the nuclear force in this framework [!, 0, H, @] but some questions regarding the 
power counting remain 0, H, @, E3, II 



The quark-mass dependence of the chiral nucleon-nucleon (NN) interaction was studied 
in Refs. [H, [H, [l4|. These studies found that the inverse scattering lengths in the Si— 



3 Di and 1 So channels may both vanish if one extrapolates away from the physical limit to 
slightly larger quark masses. 1 Subsequently, it was pointed out that QCD is close to the 
critical trajectory for an infrared RG limit cycle in the 3-nucleon sector. This led to the 
conjecture that QCD could be tuned to lie precisely on the critical trajectory through small 



changes in the up and down quark masses away from their physical values [15 . 

In the vicinity of this critical trajectory another EFT becomes useful. In this EFT 
observables are calculated as an expansion in powers of R/\a\, where R is the range of the 
two-body potential, and a is the two-particle scattering length. In nuclear physics this is 



the "pionless" EFT; EFT(/) [3j, M, U3, UM, llE l2Q| . EFT(/j has been used to compute a 
number of NN system observables as a function of R/\a\ (3, 0]- In this theory nucleons 
are described as point particles with zero-range interactions whose strengths are adjusted 
to reproduce the scattering lengths at and a s . The effective ranges r s , r t and higher-order 
terms in the low-energy expansions of the phase shifts are treated as perturbations. Since 
the spin-singlet and spin-triplet np scattering lengths of a s = —23.8 fm and a t = 5.4 fm 
are significantly larger than the iViV-interaction's range R ~ 2 fm, EFT(/)'s expansion in 
R/\a\ is a useful tool for analyzing the NN system at energies <C 1/(MR 2 ). In fact, the 
utility of this EFT is not confined to nuclear physics, since systems where R/\a\ is a small 
parameter also occur in atomic, molecular, and particle physics. In the limit that R/\a\ — > 



1 Due to the nuclear tensor force, the 3 Si and 3 Di channels are coupled. This mixing is included in the 
%EFT calculations while it appears as a higher-order effect in the pionless EFT discussed below. For 
simplicity, we will only refer to the 3 Si and partial waves in the following. 
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such systems exhibit "universal" features, which are completely independent of details of 
the inter-particle potential [2l|. The calculations of Refs. [H, [H, ll| suggested that NN 



systems with a common light-quark mass m u = rrid chosen to yield m n ~ 200 MeV would 
be close to this universality limit. 

A particularly striking feature of this limit occurs in the three-nucleon system. In the 
1970s Efimov showed that if \a\ is much larger than the range R of the interaction there are 
shallow three-body bound states whose number increases logarithmically with \a\/R. In the 
'resonant limit' a — > ±00, there are infinitely many shallow three-body bound states with 
an accumulation point at the three-body scattering threshold. If the particles are identical 
bosons, the ratio of the binding energies of successive states rapidly approaches the universal 
constant Aq ~ 515. Efimov also showed that low-energy three-body observables for different 
values of a are related by a discrete scaling transformation in which a — > X^a, where n is an 



integer, and lengths and energies are scaled by the appropriate powers of A |22|, [23J . The 
discrete scaling symmetry is the hallmark of an RG limit cycle [24| . 

The Efimov effect can also occur for fermions with at least three distinct spin or isospin 
states and therefore applies to nucleons as well. The extension of EFT(^)to the three- 
nucleon system allows a systematization of Efimov's "qualitative" approach to the three- 
nucleon problem (25[. In EFT(^), treatment of the S1/2 nd partial-wave requires the inclu- 
sion of a three-body force at leading order in the power counting j26[. It is this three-body 
force (commonly denoted Hq), whose renormalization-group evolution is governed by a limit 
cycle. Thus, once one piece of three-body data is used to fix the three-body force at a 
given regularization scale, other three-body observables can be predicted in EFT(/). This 

provides an explanation for correlations observed empirically in the three-nucleon system, 
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e.g. the "Phillips line" correlation between B 3 and the doublet nd scattering length a nd . 

But the results such as the Phillips line which were discussed in Ref. |26| were only 
obtained there up to corrections of order R/\a\ and kR, where k is the typical momentum 
of the process under consideration. The corrections to three-nucleon scattering observables 
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linear in R/\a\ were considered in Refs. |27|, [28|. Efforts to compute (R/\a\) and (kR) 



corrections to three-nucleon observables were pursued in Ref. 23|. However, recent work 
within a reformulation of the equations describing three-nucleon scattering shows that, once 
one three-body datum is used to fix the leading-order three-nucleon force, the inclusion of 
corrections up to and including (-R/|q|) 2 effects is straightforward and does not require any 
additional three-body input 30|, |3l|. (For a general analysis of the order in R/\a\ at which 
three-body input is first needed in a given nd partial wave, see Ref. [3~2|.) 

These efforts open the way for precision calculations of three-body observables in 
EFT(/) . Calculations of three-body observables in this EFT are much simpler than in 
the xEFT, and the computational effort is significantly smaller. This is particularly so near 
the critical trajectory, where the xEFT is trying to bridge a huge range of distance scales: 
from a short-distance scale 1/A X of fractions of a fermi to an NN scattering length of lit- 
erally hundreds of fermis. The pionless EFT provides a natural way to build in the physics 
that occurs between the distance scale R and the unnaturally large scattering lengths that 
occur in this regime. 

Therefore the xEFT and EFT(^) offer mutually complementary approaches to an under- 
standing of the quark-mass dependence of few-nucleon system observables. In order to use 
EFT(/) to understand quark-mass dependence we, of course, need input from the %EFT, 
but once that input is in hand predictions for a variety of few-nucleon-system observables 
can be derived up to a given order in the R/\a\ expansion. 
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A first exploratory study along these lines was carried out in in EFT(^) in Ref. 15] . The 
quark-mass dependence of the nucleon-nucleon scattering lengths from Ref. [3] was used 
as input to that calculation. In Ref. 33j, a detailed investigation of the possibility that 
1/a could vanish in both the 3 Si and 1 S channels, at the same was performed using an 
NN potential computed up to next-to-leading order (NLO) in the chiral expansion. The 



"critical" pion mass studied most thoroughly in Ref. |33j was nf™ 1 = 197.8577 MeV. For that 
case the NLO ;\EFT calculation was matched to a leading-order (LO) EFT(^) computation. 
The m^-dependence of 1/a in both channels, as well as the pion-mass dependence of the 
triton binding energy, was used as input to that computation, and the spectrum of excited 
three-body states for m, close to rn c ^ %t was computed. At LO that spectrum consists of an 
infinite tower of states with binding energies that obey: 

" :i ~ 515.035, (1) 



d("+1) 
n 3 



where B$ is the binding energy of the nth excited state in the tower. This LO prediction 
will be modified by higher-order terms in EFT(^) which are larger for the more deeply 



bound states. Examining the 'data' of Ref. [33j for this ratio we see that it is 543 for n = 



and 5.2 x 10 2 for n = 1. Here we show that EFT(/) computations of the ratio B% / B^ 
converge rapidly for all n ^ 0, with the prediction (jTJ receiving corrections of at most 0.1% 
at NLO, and vanishingly small corrections at orders beyond that. 



This extends the work of Ref. [33j, because we use the same NLO ;\EFT calculation 



as input, but compute observables up to next-to-next-to-leading order (N 2 LO) in EFT(^). 



The work of Refs. [30j, [3JJ shows that in order to do this the only additional information 
we need is the effective range in the 3 Si and 1 So channels as a function of m n . With these 
two effective ranges, together with knowledge of l/a s (m^) and Bjjn^) and 5 3 (m^), we can 
use EFT(/) to predict the behavior of the excited states of the three-body system and the 
doublet nd scattering length ajj as functions of m n in a domain of around m c ™ 1 . We 
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also provide NLO predictions for the quartet nd scattering length a nd as a function of m^. 

In the process of providing these predictions we are presented with the opportunity to 
examine the convergence pattern of EFT(/) in the three-nucleon sector. Since consideration 
of different pion masses yields different values of R/\a\, EFT(^) will not converge at the same 
rate at all values of the pion mass. Furthermore, the presence of multiple three-nucleon 
bound states near the QCD critical trajectory provides us with the opportunity to examine 
how this EFT converges for bound states with different binding energies. In fact, we will find 
that near the critical trajectory the limiting factor in the accuracy of the calculations comes 
from corrections which scale as kR (with k the momentum characteristic of the bound state 
being studied) and not from the often-quoted expansion parameter R/\a\, which is actually 
very small for near m c J lt . 

This paper is structured as follows. In Section [Til we review the pertinent results obtained 
in chiral effective field theories for quark-mass dependence of observables in systems with 
A — 0, 1, 2, and 3. In Sec. IIII[ we explain how these results are used as input to equations 
that describe neutron-deuteron scattering (and neutron-deuteron bound states) up to N 2 LO 
in EFT(/) . We present the results of our N 2 LO computation in Sec. IIVI together with a 
discussion of their convergence pattern and conclude in Sec. |V] 
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II. CHIRAL EFFECTIVE FIELD THEORY 



In order to understand the quark-mass dependence of hadronic observables we need a 
theory that encodes the soft breaking of QCD's SU{2) L x SU{2) R chiral symmetry by the 
quark mass term in the QCD Lagrangian. Chiral effective field theory is a low-energy 
theory with the same (low-energy) symmetries and pattern of symmetry breaking as QCD. 
As such it provides a systematic way to compute the quark-mass dependence of observables 
in few-nucleon systems. 

In the A = and A = 1 sectors the technology to do this is well-established, coming 
under the name "chiral perturbation theory" 0]. Scattering amplitudes are expanded in 
powers of the small parameter: 

P= m ^ P (2) 
Air f n , m p 

where p is the typical momentum of the process under consideration. Since, according to 
the Gell-Mann-Oakes-Renner relation, m 2 ~ m q , chiral perturbation theory provides access 
to quark-mass dependence, up to an accuracy that is determined by the order to which the 
computation is carried out. 

In the A = sector many observables have now been computed up to two loops, which 
is equivalent to a computation up to effects of 0{m^) = 0(raj?). In the case of A = 1 most 
computations have been performed up to "complete one-loop order" , which means that they 
include all effects up to 0{m\) = 0{m 2 q ). 

For A > 2 the presence of non-perturbative effects makes the chiral counting more inter- 
esting. Weinberg proposed that the chiral perturbation theory expansion could be applied 
to the nucleon-nucleon potential in the A = 2 system, and, more generally, to the over- 
all nuclear potential in a many-body system In this counting three-nucleon forces do 
not occur until next-to-next-to- leading order in the chiral expansion, and are suppressed by 



three powers of the small parameter [34j]. The calculations of Refs. [15|,j33j were based on a 



nuclear force computed up to O(m^) = 0(m q ) in the chiral expansion [35j 

In this section, we review the relevant results on the quark-mass dependence of hadronic 
observables that entered that calculation. Since the results given there for the NN and NNN 
systems are only valid up to effects of 0(m q ) we in general do not provide expressions that go 
beyond this accuracy. We will see that present uncertainties in the quark-mass dependence 
of observables in the A = 2 system preclude reliable predictions for the m g -dependence of 
three-body binding energies and scattering lengths. However, as long as there is a domain 
of m q values where both NN scattering lengths become large with respect to then the 

approach of the Section UTT| that uses EFT(/), together with a small amount of input from 
xEFT, to predict the quark-mass dependence of a variety of A = 3 observables, will always 
be applicable. This is true irrespective of the precise value of m^™* — or even whether the 
critical trajectory is only approximately realized. Therefore, EFT(/) can always be used to 
propagate information from the xEFT through to few-nucleon system observables — as long 
as \a\m v 3> 1 and p/m^ <C 1. And, as the accuracy of the ;\EFT input to EFT(/) improves, 
the results obtained using the expansion around the universal limit \a\ = oo will become a 
better approximation to the quark- mass dependence of full QCD. 

a. A = At the order we work to here, quark-mass dependence is synonymous with 
pion-mass dependence because of the Gell-Mann-Oakes-Renner relation: 

ml = -(m u + m d )(0\uu\0)/fi, (3) 
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where (0|mt|0) ~ (—225 MeV) 3 is the quark condensate, and f n = 92.4 MeV is the pion- 
decay constant. In the following, we will therefore discuss all our results in terms of pion-mass 
dependence. This is more convenient for nuclear applications, since the pion mass is easy to 
adjust. The corresponding change in quark masses can then be read off from Eq. ([3]). 

b. A = 1 The nucleon mass' dependence on has now been calculated to complete 
two-loop order [36j. However, here we need only the result up to 0(m q ). This can be written 
as: 

M = Mo-4c im l + 0(ml), (4) 

where the leading non-analytic piece at 0(m\) is easily calculated from the one- loop nucleon 
self-energy (in dimensional regularization with minimal subtraction), and the LEC C\ is 
related to the nucleon u-term. Meanwhile, Mq is the nucleon mass in the chiral limit. If we 



take ci = -0.81 GeV -1 |37[, we have M = 880 MeV. Correspondingly, M(m n ) increases 
to over a GeV at the putative critical pion mass m n ~ 200 MeV. Epelbaum et al. claim 
that this has only a small effect on spectra and so do not consider the effect of this change 
of m G3, Hi . Below we will provide an independent assessment of how this change in M 
affects the results near = m c ^ lt . 

The other A = 1 quantity of relevance for our study of nuclear bound states is the 
pion-nucleon-nucleon coupling. For this we have the result: 

9 * NN ^fl + 2A-^ 18 V (5) 



M U\ g A 

where qa = 1.26 is the physical value of the nucleon's axial coupling. This is the Goldberger- 
Treiman relation in the chiral limit, and here it is supplemented by an 0(m q ) term involving 
the LEC dja. Using a value of g n NN extracted from a phase-shift analysis of low-energy 7rN 



data [38|, |39J we get di$ = —0.97 GeV~ [l4|. Also appearing in Eq. (jSJ) is the ratio A, which 
encodes the m g -dependence of the ratio gA/ fir- A is just the fractional change in this ratio 
at an arbitrary value of compared to the value at the physical m n , m^ hys , i.e.: 

A = (9A/fn)m n ~ (g A / fn) m Phys 
{9 A / fn) m Phv* 

The expression for A in terms of low-energy constants appearing in and C^l can be 



found in Ref. 33 



c. A = 2 The quark-mass dependence of the chiral NN interaction was calculated to 



next-to-leading order (NLO) in the chiral counting in Refs. [12|, [13|, Il4j . In addition to the 
quark-mass dependence of the pion mass and nNN coupling that appear in the one-pion- 
exchange potential the short- distance part of the force acquires a quark-mass dependence 
of its own. This is essential for correct renormalization [40] . This quark-mass dependence 
can be described by parameters D s and D t , so that the NN potential in a particular S-wave 
(t=3 Sl , s =%) is: 

V = C s/t + D s/t [ml - (mf ys f] + V OPE + V<$ E , (7) 

where Vope is the (partial- wave projected) one-pion exchange with the already discussed 
variation in g^NN dp a nd ([3]) taken into account, and V^p E is the well-known "leading" 



two-pion exchange [35j, |4l|, |42| . Note that to the order considered here it is only necessary to 



take into account the m„--dependence of g n NN (or, equivalently, gA and /„.) in the LO piece 
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f 2 A 2 D 

Oi s t = — — (8J 



of the NN potential. We do not take into account the m n - dependence of qa and f w in the 
two-pion-exchange potential but instead simply use their physical values there. Considering 

(2) 

only "explicit" dependence in Vj.p E is sufficient at the order we work to here. 

Meanwhile, in Ref. [33] the unknown parameters D s / t were converted into dimensionless 
numbers via: 

x2 A 2 T, 

>s/t 

16n 

At present next-to- nothing is known about the values of a. But predictions for the pion- 
mass dependence of S-wave scattering can be obtained by varying the a's within naturalness 
bounds. The uncertainty in the behavior of A as a function of m n then further enlarges the 
range of results obtained for, say, the A^ A" scattering lengths as a function of m n . 

Because of these uncertainties, the study of Ref. [33] did not attempt to predict the 
precise value of the critical pion mass where the inverse scattering lengths in both S-waves 
vanish. This strategy is motivated not just by ignorance as regards the two LECs appearing 
in the a's, but also because we have not considered the impact that m u ^ rrid would have 
on the A^A" potential. Such effects must be included if the Braaten-Hammer conjecture, 



as originally formulated, is to be investigated [15J. Instead, in Ref. [33] a parameters were 



chosen in such a way that a critical trajectory occurred. For critical pion masses in the 
range 175 MeV < m^"* < 205 MeV, it is always possible to find an appropriate set of a 
parameters. Three sets of a parameters were given in Ref. [33[ , with each one corresponding 
to a different critical pion mass m c ^ lt . The computations of Ref. [33[ focused particularly on 
a choice at = —2.5, which yields a s = 2.138598 and a critical pion mass m^ lt = 197.8577 
MeV. Although it is unlikely that physical QCD will correspond to this particular solution, 
many aspects of the limit cycle are universal and do not depend on the exact parameter 
values [2l|]. Therefore in what follows we take the results reported in Ref. [33] for this 
choice of xEFT parameters as both the input data for EFT(/f) and the experimental data 
to which we will compare our EFT(^) results for the pion-mass dependence of three-nucleon 
observables in the vicinity of the QCD critical trajectory. 

As already discussed in the introduction, the A^/V inputs we need from the xEFT for our 
NNLO EFT(^) calculation are scattering lengths and effective ranges in the spin-triplet and 
spin-singlet channels. The scattering lengths at and a s can, of course, be exchanged for the 
momenta jt and 7 S characterizing the position of the bound-state/virtual-state poles in the 
two-body propagator. If higher-order terms in the effective-range expansion are neglected, 
the relation between the two quantities is simply 



- (l - yfl - 2r a /a a ) , (9) 



7a 

with a — s,t indicating a particular A^A^ channel. The advantage of using 7 a , rather than a a , 
is that the A^A^ pole position is not modified by higher-order corrections. In the following, 
we will therefore use 7 t (m 7r ), 7< ! (m 7r ), r^m,,-), and r s (m n ) as the AW input. 

The behavior of the input pole momenta 7 s (m 7r ) and 'jti'^n) i n the vicinity of m n = m c ™ t 
is depicted in the inset in the left panel of Fig. [TJ As promised, the pole momenta vanish 
at the critical value of the pion mass rrf™ 1 = 197.8577 MeV. For pion masses below m^ lt 
74(777.7,-) > and the deuteron is bound. As 74(777^) decreases, the deuteron becomes more and 
more shallow and finally becomes unbound at the critical pion mass. Above the critical pion 
mass the deuteron exists as a shallow virtual state. In the spin-singlet channel, the situation 
is reversed: the "spin-singlet deuteron" is a virtual state below the critical pion mass and 
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FIG. 1: Left panel: Pole momenta in the 3 Si (dashed line) and 1 So (solid line) nucleon-nucleon 
channels as a function of the pion mass (from Ref. 33|). The vertical dotted line indicates the 
critical pion mass m£™* while the inset shows the critical region in more detail. Right panel: same 
for the effective ranges (from Ref. j^). 



becomes bound above. In the right panel of Fig. [TJ we show the corresponding dependence 
of the spin-singlet and spin-triplet effective ranges j43|. The pion-mass dependence of the 
effective ranges is weak and both are close to 2 fm throughout the region displayed. In 
the region of interest r s (m n ) decreases monotonically as the pion mass is increased whereas 
r t (m v ) increases. At a pion mass of about 202 MeV the curves cross. 

which are presented in Fig.[TJdo not 



Above we pointed out that the results of Refs. [3J 
include effects due to the pion-mass dependence of the nucleon mass. Such effects can easily 
be accounted for though, by noting that the computations presented in Fig. [TJ are done in a 
non-relativistic framework, and no 1/M corrections to the AW potential are incorporated. 
Consequently the nucleon mass can be scaled out of the problem by an appropriate choice 
of units. The critical pion mass is then written as: 



rn 



crit 



0.2107M, 



(10) 



where we have employed an (isospin-averaged) nucleon mass of M = 938.9 MeV. A value 
of rn c ^ lt in MeV that includes the impact of M(m n ) can be obtained from Eq. ffTUl) by 
inserting the expression (jl]) for M(m^"*) in place of M on the right-hand side, and solving 
the resulting equation for rrf™ 1 . This yields rrf™ 1 = 220 MeV, with the accuracy limited to 
two digits by the precision with which c\ is known. 

We can use the same argument to assess the impact of the m^-dependence of M on two- 
body observables when ^ mf"*. Incorporating the pion-mass dependence of the nucleon 
mass in the calculation of Refs. [33|, |43[ in this way would thus lead to a change of the scale 
on both axes in the left panel of Fig. CD by a factor of M(m 7r )/M. The x-axis in the right 
panel would be rescaled by the same factor, while the y-axis there would be rescaled by 
M/M(m n ). The effect of such rescaling is largest at the highest pion masses considered, 
but even there it is less than 15%. Other uncertainties in the calculation of the functions 
7s("i.7r), lt( m n), T' s (m 7r ), and r t (m^) are larger than this, and so from here on we follow 
the procedure of Refs. [i 
dependence of M. 



43] and ignore effects on these quantities due to the pion-mass 



S 



d. A =3 We also need one three-body datum as input to EFTJY) in order to perform 
our N 2 LO computation of three-body system observables. In Ref. 33|, the binding energies 



of the triton and the first two excited states in the vicinity of the limit cycle were obtained 
from the solution of the Faddeev equations for the NLO ;\EFT potential. For much of the 
pion-mass range of interest to us there is only one bound state, and so, in that region, 
100 < < 190 MeV, we renormalize EFT(^) using the binding energy of the triton 

ground state, B$ , and then predict the nd doublet scattering length, ajj- But our results 
for pion-mass dependence are particularly interesting in the "critical region", which we 
define to be values of for which the triton has at least one excited state. For the choice 
m crit _ 197 8577 MeV the calculations of Ref. 33( indicate a critical region corresponding 
to 190 MeV < m n < 210 MeV. The first excited state that is present in this region is 
much shallower than the ground state, and is transparently within the domain of validity of 
EFT(/) . Thus errors in that EFT are smaller if we renormalize to its binding energy, Btp, 
rather than to Bf\ The EFT(^) result for ^ 0) then becomes a prediction and we can also 
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predict a n ' d . The ;\EFT results for both the ground and first-excited state in the "critical 
region" are shown in Fig. [2], together with our N 2 LO EFT(/) calculation. 



III. THE PIONLESS EFT TO N 2 LO 

In the critical region where there is more than one three-nucleon bound state, the length 
scales in the three-nucleon problem change rapidly as is varied. The binding energy of 
the deepest lying three-nucleon state changes slowly, and is still within a factor of two of 
its experimental value when m c ™ 1 is reached. However, in this region the neutron-deuteron 
scattering length varies between — oo and oo, and does so an infinite number of times. 
A single cycle between -oo and oo takes place over the range of pion-mass increase that is 
needed to make an additional shallow three-nucleon ("Efimov") state appear as — > m c £ lt . 
Consequently, the numerical effort which has to be devoted to the calculation of three- 
nucleon observables increases significantly once is large enough that the first Efimov 
state appears — or is close to appearing. We advocate determining the dependence of key 
two- and three-nucleon observables on the pion mass with %EFT, then using those results 
as input to the simpler EFT(/) and employing that EFT to analyze the behavior of other 
three-nucleon observables in the critical region. 

EFT(/) is a systematic expansion in contact interactions where — for momenta k ~ 1/a — 
the small expansion parameter is R/a. The corresponding low-energy Lagrangian for the 
neutron-deuteron system is given by [26| 



+g t (t\N T T 2 (Ji(J2N + h.c) + g s [s)N T a 2 a j T 2 N + h.c^ 
-G 3 N* (glitiarftvav + ^^[(i^)^- + h.c] + g^SjT^ s f r^j N + ... , (11) 

where N represents the nucleon field and U(sj) are the di-nucleon fields for the 3 Si ( 1 So) 
channels with the corresponding quantum numbers, respectively. The Oi (tj) are Pauli 
matrices in spin (isospin) space, respectively, and the dots indicate additional terms with 
more fields/derivatives. 
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To renormalize three- nucleoli observables in the nd S1/2 channel jvithin this framework, a 



three-nucleon force symmetric under SU(4) spin-isospin rotations |26|] (here denoted by G3) 



or equivalently a subtraction has to be performed [44j. |45| . Therefore, one three-body input 
parameter is needed for the calculation of observables in this channel — the one in which the 
triton exists. 

In order to compute three-nucleon observables we need the full two-body propagator r, 
which is the result of dressing the bare di-nucleon propagator by nucleon loops to all orders. 
The EFT is arranged to reproduce the effective-range expansion for r, i.e. 

2 1 

tJE) = — ; , (12) 

V ' vrM _ 7q + vQH + if ( 7 2 + ME) V ; 

where E denotes the two-body energy in the two-body cm. frame, 7 a gives the NN pole 
position, r a is the effective range, and the index a = s, t indicates either the singlet or triplet 
NN channel. In the triplet channel the pole position jt is related to the deuteron binding 
energy = 7 4 2 /M. The form (fl~2l) cannot be directly used as input to EFT three-nucleon 
calculations because it has poles at energies outside the domain of validity of the EFT. 
Therefore, the propagator cannot be employed within a three-body integral equation which 
has cutoffs A > l/r a unless additional techniques to subtract these unphysical poles are 
implemented. Instead of using the propagator in the form above we will expand it to nth 
order in R/a: 

r M(E) - S ^( E ) d<V) 
Ta E + ji/M- {U) 

For n < 3, the residue is given by 

S£\E) = ^E(y) 1 [7. + V^ME}^ , (14) 

i=0 

while for n > 3 higher-order terms in the effective-range expansion contribute to S(E). 
The set of integral equations for the nucleon-deuteron K-matrix generated by this EFT 



(neglecting, for the moment, the nd coupling G 3 ) is given by [26l. 145 

/A o 112 

dq" q"*Z u (q, q"; E)ri n) (E - ^)K^(q", q'; E) 

/A o 111 

dq" q»*Z ts (q, q"; E)r^(E - ^)K%\q", q'; E) , 

/A o 1/2 

dq" q" 2 Z st (q, q"; E)r^(E - ^)K^\q", q'; E) 

/A o 112 

dq" q" 2 Z ss (q, if\ E)t<?\E - —^(^ 4\ E) , 



(15) 

where n is the order of the calculation (in what follows we assume n < 3) and V indicates 
a principal-value integral. The formulation in terms of the K-matrix is useful as long as we 
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FIG. 2: Triton ground and excited state energies B% in the critical region. The circles (ground 
state), squares (first excited state), and diamonds (second excited state) give the x^FT result, 
while the solid lines are calculations in the pionless theory to N 2 LO. The vertical dotted line 
indicates the critical pion mass m^" lt . The thresholds for three- nucleon states to be stable against 
breakup into a single nucleon plus an AW bound state are given by the dashed lines. The zero of 
energy corresponds to the opening of the 3N channel. 



are only interested in bound states or nd scattering below the three-nucleon threshold. The 
Born amplitude Z a p is given by 



Z a p(q, q;E) = -\ a p— log 
qq' 



ME\ 

q 2 - qq' + q' 2 - ~ME J 



(16) 



with the isospin matrix 



A, 



(17) 



1 

\a={s,t};f3={s,t} - J 

Below three-nucleon breakup threshold the K-matrix is related to nd phaseshifts through 



1 -3 
-3 1 



K(k,k;E) 



3M 
8^k 



tan 5 . 



(18) 



In Refs. HSU, it was shown that the use of one subtraction allows rW up to N 2 LO (n = 2) 
to be inserted in Eqs. ( fT5j) with cutoff-independent predictions resulting. This is in marked 
contrast to the unsubtracted equations written above, which generate cutoff variation of 



461 ] . The price that is paid for 



0(1) in low-energy observables if left unrenormalized [2i 

1/2 

the improved ultraviolet behavior of the subtracted equations is that a nd appears as an 
input parameter in the integral equation. For further details on this subtraction method, we 
refer to Refs. El EH 



In the next section we will use the subtracted equations to perform 
calculations in the pionless EFT up to N 2 LO. 

To close this section we note that the Lagrangian (jlip is symmetric under Wigner SU(4) 
spin-isospin rotations provided g s = g t and A s = A t . Moreover, the integral equations 
( TT5|) have an approximate SU(4:) symmetry for momenta q 3> 7 s ,7t, even if g s ^ g t or 
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FIG. 3: Binding energy Bj, of an excited state of the triton as a function of 7 S and in the vicinity 
of m£ . The computation is performed here to leading order in EFT(/f) . 



A s ^ A t . (For a discussion of this symmetry for the NN system see Ref. (47]].) Because 
of this symmetry, an S , f/(4)-symmetric three-body force is sufficient for renormalization 
26). In Fig. El we illustrate a manifestation of this symmetry in the bound-state spectrum 



calculated at LO in EFT(^). We show the binding energy of a triton excited state in the 
critical region as a function of 7 S and jt- In the lower left corner, the considered triton state 
does not exist. As 7 S and jt approach zero, the state appears at threshold and then becomes 
more and more bound. The figure is clearly symmetric under reflection on the main diagonal 
corresponding to a specific spin-isospin rotation that exchanges 7 S and 7 t . 

Beyond leading order in EFT(^) S'f/(4)-breaking effects enter the computation through 
the difference in singlet and triplet effective ranges r s — r t . The right-hand panel of Fig. [1] 
displays the interesting feature that — at least for the choice of %EFT short-distance coeffi- 
cients being studied here — r s — r t vanishes at a pion mass only a little larger than m ™ 1 . As 



remarked upon in Ref. 33] , this results in SU (4)-breaking effects near the critical point that 
are much smaller than one might naively expect, since r s — r t <C R, the range of the NN 
interaction. More generally, Fig. [1] suggests that a suitable choice of a s and a t could lead to 
r s — r t = 7 S = 7f = at a single U SU (4) critical pion mass" . It is unclear whether a xEFT 
that has this feature has anything to do with QCD. But the EFT(/) valid near a critical 
trajectory which also has r s = r t could be built using S'L r (4)-symmetric operators for the 
R/a expansion, with S'6 r (4) breaking additionally suppressed by another small parameter. 



IV. RESULTS 

In this section, we solve the once-subtracted version of the integral equations (1151) for pion 
masses in the range 100 MeV < m n < 200 MeV. Two-nucleon pole positions and effective 
ranges for the relevant xEFT solution were displayed in Figfl], and in the critical region we 
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FIG. 4: Binding energy of the triton ground state in the critical region computed in the pionless 
EFT at LO (dashed), NLO (dash-dotted) and N 2 LO (solid). The dotted line gives the (interpo- 
lated) xEFT results for comparison. 

take as our one three-body input the energy of the first triton excited state B^. We provide 
a detailed discussion of the critical region, and also examine the behavior of EFT(^) for pion 
masses between the critical region and the physical value. In the latter case the three-body 
datum chosen as input is the triton binding energy B% . 



A. Bound-State Spectrum 

Turning first to the critical region, we display the results obtained from solving the 
homogeneous version of Eq. f|T5|) for negative energies. The spectrum of triton states as a 
function of m n at N 2 LO is shown in Fig. [21 and compared to the xEFT result. 

The binding energies from xEFT are given in Fig. [2] by the circles (B^), squares ), 
and diamonds (B^) {33|. The dashed lines indicate the neutron-deuteron (m n < m c J lt ) 
and neutron-spin-singlet-deuteron (m n > m c J' lt ) thresholds. For B3 less than these energies 
three-body states are unstable against breakup into the 1 + 2 configuration. Directly at the 
critical pion mass, these thresholds coincide with the three-body threshold and the triton 
has infinitely many excited states. The solid lines show our N 2 LO calculation in the pionless 
EFT where the first triton excited state was used as input. Our results for the second excited 
state and the ground state reproduce the xEFT results very well. An important point here 
is that the binding energy of the triton ground state varies only weakly over the whole range 
of pion masses. Indeed £?g (m£T**) ~ 0.5B^\m^ hys ). The excited states are influenced by 
the 1 + 2 threshold and their energies vary more strongly. 

In order to study the convergence of the EFT(^) more thoroughly we now look at the 
prediction for the triton ground-state energy in some detail. In Fig. HJ we show our results 
for the binding energy of the three-nucleon ground state for pion masses in the "critical 
region" The dashed, dash-dotted, and solid lines denote, respectively, the LO, NLO, and 
N 2 LO results for the triton ground-state energy, while the dotted line gives the (interpolated) 
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(2) 



FIG. 5: Diagrams contributing to the NLO shift in the binding energies. The shaded square 
in Diagram 1 represents a non-trivial effective-range insertion. The shaded circle in Diagram 2 
denotes the insertion of the subleading three-body force H\. The shaded blobs are LO nd vertex 
functions. 

xEFT result of Ref. [33|]. The leading-order EFT(^) results describe xEFT's Bf'im^) curve 
surprisingly well, although the shape of the m^-dependence is clearly wrong. The NLO 
computation does not seem to yield significant improvement in the overall agreement. The 
results of our N 2 LO calculation are on average better than the LO results. 

At first sight this slow convergence is surprising, since the parameter — > at = 
rrf™ 1 . But this simply means that the slow convergence in the critical region must be due 
to the significant (compared to 1/r) binding momentum k = ^/4MB 3 /3, which is about 65 
MeV in this pion-mass range. 



B. Understanding Range Corrections in EFT(/f) 

A rough estimate of the impact of higher-order corrections on EFT(/) predictions can be 
obtained by applying naive dimensional analysis. The relevant scales for the orders under 
consideration are 7, r and the binding momentum k. (For simplicity we will drop the indices 
s and t in the ensuing discussion since the scaling with r s and r t is the same.) It follows 
therefore that the error for an arbitrary three-body observable O3 at order j should scale as 

— -A j+1 (^-Y +1 + B J+1 M^, (19) 



3 ^ l \2, 

where A,+i and Bj + \ are (observable-dependent) numbers of order one, and k is the mo- 
mentum associated with 3 . 

However, this naive estimate clearly fails to explain the pattern of convergence displayed 
in Fig. [2j A more thorough understanding of this pattern of corrections can be gained 
by examining the way that the NLO correction at the critical pion mass scales with the 
parameters k and r. We will analyze a perturbative expansion of EFT(^) about the LO 
result (281]. When m % = m^™*, 7 = and 7r corrections do not exist. In this "scaling limit", 
the two-body propagator at NLO simplifies to 



7rM 



1 1 

+ 



(20) 



In Fig. [5] we display the diagrams which have to be taken into account for analysis of nr 
corrections. Diagram (1) represents a perturbative insertion of r/2, a la the second term 
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of Eq. (1201) . Diagram (2) denotes the inclusion of the subleading three-body force Hi. The 
presence of H\ is necessary in order to absorb a divergence produced by diagram (1), but 
H\ always appears in a fixed combination with Hq, i.e. it does not need to be determined 
by any additional three-body data |2a |. 

An analysis of the superficial degree of divergence of these diagrams gives amplitudes 
which scale with momentum p as: 

2M' M 2 K ' 

However, the only remaining momentum scale in the problem is the binding momentum for 
the bound state under consideration, i.e. K^ n \ The NLO shift in the binding energy of the 
nth bound state is therefore 

ABi n) = z(a4-^-(3^^ 2 ) , (22) 
3 V 2M ' M 2 ) ' y ' 

where a and (3 are numbers of order one, which, due to the discrete scale invariance of the LO 
spectrum in EFT(^) , are independent of n. Z denotes the wave-function renormalization, 

and a similar diagrammatic analysis leads us to conclude that Z ~ . We therefore 

obtain 

Afl, W =(^ («•"»)• («<•>)') ■ (23) 

If we have renormalized to the experimental binding energy of the mth bound state at leading 
order we want to preserve that experimental value at NLO, and so we demand AB^ m) = 0. 
This fixes the strength of Hi to be 

arM 1 

2(3 ^ 

and the shift in the binding energy of the nth state therefore becomes 



ft = • < 24 > 



^T = 6 — i 1 "^))- < 25 > 

Parametrically, the fractional change in is of order /c^ n V, as expected. However, 
the number in brackets in Eq. (1251) is not necessarily of order 1, since the ratios of binding 
momenta are large. Because of the discrete scaling symmetry of the three-nucleon spectrum 
which is realized in our LO EFT(^) calculation we have 

= (22.7)—. (26) 

If the state which has been used as a renormalization point (m) is more weakly bound than 
the states for which we are making a prediction (n) we have m > n, and so the NLO shift 
will be amplified by a factor of (22.7) m_n . This to some extent explains the rather large 
shift in B^ when we renormalize to B% . Conversely, if the state used for renormalization 
is more strongly bound then the state under consideration, and so m < n, the NLO shift 
will be 

AB^ K Wr , x 

— ss a——, 27 
B {n) 2 y ' 
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[MeV] 




Bf ] [MeV] 




LO 
NLO 
N 2 LO 


3.579 
3.369 
3.594 


515 
485 
517 


1.349 • 10~ 5 

1.350 • 1(T 5 
1.350 • 10~ 5 


515.0 
514.8 
514.8 


X EFT 


3.774 


543 


1.329 • 10" 5 


523.1 



TABLE I: Binding energies B^ and B^ and ratios B^ /By and B^ / 'B^ at the critical pion 
mass. The three-body input parameter has been adjusted such that the xFFT result for the first 
excited triton state B^p = 0.00695 MeV (at m n = m^~ lt ) is reproduced. The first three rows show 
the results in the pionless EFT at LO, NLO, and N 2 LO, while the last row gives the xFFT result 
from Ref. 33]. 



in accordance with naive dimensional analysis. 

Equation ( 1231) predicts that EFT(/) should converge more smoothly for states which 
are less bound than the state which has been used for renormalization. So we now turn 
our attention to the convergence of predictions for binding energies of excited states of the 
triton, where we expect EFT(^)to work much better than it does for the ground state. 

Since we used as input, the simplest observable we can test this hypothesis on is the 
ratio BP /By . (Later we will examine the doublet nd scattering length, a^ 2 , but there we 
do not have ;\EFT data with which to compare.) Using Eq. (123|) . and remembering that we 
are effectively choosing Hi such that ABP = , we expect 



BP ( k& 



J 3 

o(2) 
n 3 



>r 



515 x 1-5— - + ... , (28) 



with the binding momentum of the second excited state. Since < 0.1%, this 

suggests that By / By will change by less than 0.1% from LO to N 2 LO. 

In Table H] we show the ratios of the first three Efimov states extracted from the xEFT and 
pionless EFT calculations. In Ref. [33[ values of 542.9 and 523.1 were given for B^/BP 
and BP/By respectively. Given Eq. ( 1251) we expect that EFT(^) converges slowly for 
binding momenta on the order of k^°\ so we do not anticipate reproduction of By /BP 
with accuracy better than a few per cent at N 2 LO. This expectation is borne out by the 
results of Table HI 

In contrast, ranee corrections are miniscule in the case of B ( p/Bf\ We have calculated 
NLO and N 2 LO kt corrections to this ratio at m, = m c ^ lt . We find a shift at NLO in the 
ratio that is numerically significant, and on the order of the expected 0.1%. No N 2 LO shift 
in bP /BP is seen within our numerical accuracy. We summarize this result as: 

B {i) 

-|- = 514.8 ± 0.02 ± 0.0004, (29) 
B? 

where the first error is from the numerical precision of our calculation, and the second is from 
higher-order terms in the EFT(^) expansion. We note that the xEFT calculation of this 
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FIG. 6: Doublet neutron-deuteron scattering length a ' d in the critical region computed in the 
pionless EFT. The solid line gives the LO result, while the triangles and circles show the NLO and 
N 2 LO results. The dotted lines indicate the pion masses at which aj^ diverges. 



quantity is numerically quite delicate, since it requires a momentum mesh that bridges scales 
ranging from A = 540 MeV (the cutoff in the chiral potential) to the binding momentum 
« 130 keV — more than three orders-of-magnitude difference. Indeed, reexamination of 
the calculation of the second excited-state energy in Ref. [33[ showed that the numerical 



value of 523.1 quoted there is only accurate to two significant digits [481 ]. There is thus no 



disagreement between our result and that quoted in Ref. [33]. This emphasizes the 
ability of EFT(^) to easily obtain high-precision results for three-body observables. 

Therefore, at the critical pion mass, Bjp /B$ differs from the Efimov prediction ([T]) by 
less than 0.1%. This allows us to predict that the higher excited states will obey 



B 



(n+l), 



rn 



crit\ 



(515) n . 



(30) 



for all n > 1. This result, which is a rigorous consequence of the discrete scaling symmetry 
of the Efimov spectrum and the scaling of the corrections to that spectrum, could be a useful 
constraint on the infinite tower of excited states if the critical trajectory is ever realized in 
a lattice QCD computation. 



C. Doublet Scattering Length 

Scattering observables have not been calculated in xEFT in the vicinity of the critical 
trajectory, but can be predicted in EFT(^) in a straightforward way. In principle the sub- 
tracted version of Eq. (ITB1) can be solved at any positive energy to obtain nd phaseshifts. 
In practice additional numerical techniques are needed to do this above the three-nucleon 
breakup threshold, and at the critical pion mass this coincides with the nd threshold. How- 
ever, the technology to deal with this complication is well known. Here we examine the nd 
scattering lengths as a representative of all nd scattering observables that could be computed 
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FIG. 7: Doublet neutron-deuteron scattering length a' d around the physical pion mass at LO 
(dashed), NLO (dash-dotted) and N 2 LO (solid). Our results have been connected with splines to 
guide the eye. The dotted line gives the physical pion mass. 



in EFT(^) . The scattering length is measured at zero momentum. Since we renormalize to 
a bound state with binding momentum larger than this we do not anticipate any enhance- 
ments of the type encoded in Eq. (|25|) . If one takes into account the binding momentum of 

1/2 

the deuteron, one could argue that a n(J involves typical momenta of order 7, thus the error 
in the N fc LO EFT(;f) prediction for aj^ ( or f° r t ne quartet scattering length aJ2) should 
be, at worst, (7r) fc+1 . In Fig. El we show the doublet scattering length ajj[ in the critical 
region. The solid line gives the LO result, while the triangles and circles show the NLO 
and N 2 LO results. The dotted lines indicate the pion masses at which aJ2 diverges because 
the second and third excited states of the triton appear at the neutron-deuteron threshold. 
These singularities in o^ 2 (m r ) are a clear signature that the limit cycle is approached in the 
critical region. Fig. O shows that, as expected from our scaling arguments, the higher-order 
corrections to aj^ are very small and the NLO and N 2 LO results lie on top of the LO curve. 
Considering the specific case of the doublet scattering length at = 190 MeV, we see that 

aj^im^ = 190 MeV) = (93.18 + 0.80 + 0.14) fm . (31) 

At this pion mass we are able to match EFT(/) to B$ , and we are also fairly far away 

1 /2 

from any singular points in the function a nd (m n ). Consequently, the EFT(/) results follow 
a natural convergence pattern with the expansion parameter 77", which is ~ 0.08 at this m n . 

While the Efimov effect leads to dramatic effects around the critical pion mass, the 
pionless EFT can also be used for a similar analysis in the region between mf yB = 139 MeV 
and the lower edge of the critical region. In this pion-mass domain the triton has no excited 
states. We have therefore used the ;\EFT results for the pion-mass dependence of the ground 
(and only) three-nucleon bound state as the three-body input to our EFT(/) computation. 
We then predict the pion-mass dependence of the doublet neutron-deuteron scattering length 
aJJ. In Fig. [7]we display EFT(/) results at LO, NLO and N 2 LO. A couple of points in Fig. [7] 
are particularly worth noting. First, the convergence pattern of EFT(^) for Pd 100 MeV 
is a little peculiar, with NLO and N 2 LO corrections essentially canceling to leave the LO 
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prediction undisturbed. There is also a cancellation in the NLO correction at m n = mP hys , 

1/2 

which leaves the NLO shift in a nd there accidentally close to zero. Presumably this occurs 
because at that point the two types of NLO corrections — one proportional to the binding 
momentum and the other proportional to 77*, are equal in magnitude but opposite in 
sign. This is consistent with the observed feature that as — > m c " t the convergence 
pattern of the EFT becomes more natural, because kP't corrections begin to dominate the 
7r ones. A second interesting point is that our results suggest that a small decrease in the 
quark masses could lead to a sign change in the nd doublet scattering length. Perhaps most 
important, our N 2 LO EFT(/) computation allows us to make a firm prediction (within the 

1 /2 

scenario under consideration here): as m q increases above its "real world" value a nd will 
decrease monotonically until the first triton excited state appears, something that in the 



critical-point realization studied in Ref. [33J happens at « 190 MeV. 



D. Quartet Scattering Length 

With results for a^ 2 \m n ) in hand it is natural to ask whether EFT(^) can be used to 
say anything about nd scattering in the other channel, the quartet. In this case the integral 



equation describing the scattering is [181, |49 



K$(q, q>- E) = Z(q, q>; E)+v£ dq" q"*Z(q, q"; E)r^ (e - K$tf, q'; E) (32) 

where 

Z(q,q;E) = -—log [- — — — . 33 

2 qq' \q — qq + q — ME J 

The prefactor in Eq. ( 1331) is different from Eq. ( fT6j) . This crucial difference is a consequence 
of the spin and isospin coupling in the quartet channel, and it guarantees that the integral 
equation (1321) will yield cutoff-independent results without a subtraction being necessary. 
This is a manifestation of the Pauli principle in this channel, which precludes the appearance 
of a three-body force without derivatives. In consequence once is defined Eq. (|32|) is 
straightforward to solve. Up to N 4 LO the result for ajj[ is purely a function of j t and r t . 
Performing LO (n = 0) and NLO (n = 1) computations using Eq. fl32l) . and then expanding 
the NLO result in Taylor series about r t = 0, yields the result for this function: 

aj 2 = 1.1791 7^ + 0.5540(6) r t + O {-f t r 2 ) ■ (34) 



(For N 2 LO computations of this observable at the physical 74 and r t see Refs. |l8|, [50j.) Here 
the leading-order coefficient of 1.1791 (accurate to the number of digits written) was already 
obtained in Ref. j49|. The computation of the next-to- leading order coefficient (which has an 



error of 6 in the last digit) agrees with results from Efimov [5JJ. The result fl34|) translates 

into the results presented in Fig. E]for a^(m n ). Numerically we have, at the physical pion 

mass, 

c?J2 = 5.0911 x (1 + 0.190 ± 0.139) fm = (6.06 ± 0.71) fm, (35) 

which is consistent with the experimental value of 6.35 ± 0.02 fm for this quantity [52]. 
Indeed, the NLO calculation performed here agrees with the experimental number at the 
5% level, which is better than expected from our naive-dimensional-analysis estimate of the 
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FIG. 8: Quartet neutron-deuteron scattering length a nd around the physical pion mass at LO 
(dashed) and NLO (dot-dashed). The data points have been connected with splines to guide the 
eye. The dotted line gives the physical pion mass. The star indicates the experimental value, with 
the error bar contained within the symbol. 







size of the N 2 LO correction. The accuracy of the EFT(/) prediction for a nd {m^) should only 
improve as we move from the physical pion mass towards m c J lt . Note that once > m^ tt 

3 /2 

the quantity a nd cannot be denned, since the spin-1 NN state becomes unbound, and the 
spin-0 bound state obviously cannot be used to produce a total spin of 3/2. 



V. SUMMARY & CONCLUSION 

In this paper, we have presented a detailed study of the pion-mass dependence of three- 
nucleon observables around the critical pion mass and around the physical pion mass. We 
have performed calculations to N 2 LO in EFT(^) using the pion-mass dependence of the two- 
nucleon effective range parameters and one triton state from a NLO calculation in xEFT as 
input. 

In particular, we analyzed the convergence pattern of the effective-range corrections at 
NLO and N 2 LO in the critical region where 7 S pa j t pa 0. We found that the convergence of 
the EFT expansion is slow for the triton ground state but rapid for all remaining states in the 
three-nucleon Efimov spectrum. This behavior is expected, since the higher-order corrections 
in the critical region scale with powers of Kr, where k denotes the typical momentum scale 
of the observable under consideration. Our results demonstrate that the pionless EFT is 
well suited to describe low-energy observables in few-body systems with a large scattering 
length. For the triton ground state, it appears to be worthwhile to extend current pionless 
EFT calculations to N 3 LO, so as to further investigate the convergence pattern. For the 
triton excited states in the critical region, the N 3 LO corrections will be tiny. Although an 
N 3 LO calculation would require the inclusion of an additional three-body parameter, the 
binding energy of the three-nucleon ground state could be calculated if this parameter was 
matched to scattering data or an excited three-body state. 
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Directly at the critical pion mass, we have performed a detailed comparison with the triton 
energies from %EFT. In the zero-range limit, the three-body system displays an exact discrete 
scale invariance and the ratio of energies of two subsequent Efimov states is B$ /B^ 1 ^ = 
exp(27r/so) ~ 515.03.... For finite effective range, r, there are corrections to this number 
proportional to K^ n 'r, where «W ~ (MB^) 1 / 2 is the binding energy of the nth excited state. 
However, the exact ratio exp(27r/s ) is approached as n — > oo since «W — * in this limit. 

For the ratio of the energies of the first and second excited state, we find B$ / B^ = 514.8 
when the range corrections are included to N 2 LO. The apparent disagreement of this result 
with the direct calculation using the xEFT potential [33| was resolved by taking the reduced 



accuracy for the excited states in that calculation into account [48j. In the simpler pionless 
EFT, the excited states can be calculated to very high accuracy without much computational 
effort and the higher-order corrections are small. For the more deeply bound states on the 
other hand, the convergence of the pionless EFT is slow and high accuracy is possible in 
xEFT. This calculation provides a prime example of how xEFT and EFT(/) complement 
each other very well. 

We have also studied the pion-mass dependence of three-nucleon scattering observables. 

1/2 

Directly at the physical pion mass, the nd doublet scattering length a n ' d is unnaturally 
small. As a consequence, small variations in the pion mass lead to significant changes in 
aJ2- A s the pion mass is increased towards the critical value ajj becomes very large and 
negative, jumping to positive infinity when the first excited state of the triton appears. This 
behavior repeats as one moves closer to the critical trajectory and more and more excited 
states of the triton appear. It is a signature of the limit cycle being approached. The source 
of this strong variation of the doublet scattering length in the critical region is therefore 
understood and well described by EFT(^) . The pion-mass dependence of the nd quartet 
scattering length at the physical pion mass is much milder. As the critical pion mass is 
approached, it grows and eventually becomes infinite on the critical trajectory. In contrast 
to the doublet scattering length, this growth is simply driven by the increase in size of the 
deuteron and has nothing to do with the limit cycle. 

Our results will be useful in the context of future Lattice-QCD calculations of three- 
nucleon observables. As a first step in this direction, the NN scattering lengths have re- 
cently been extracted from NN correlators computed at quark masses corresponding to 
pion masses between 350 and 590 MeV [l|. Thanks to continuing advances in computer 
power much progress in this direction can be expected over the next couple of years. Lat- 
tice computations of NNN correlators are a high priority since they will provide access to 
three-nucleon forces directly from QCD, but algorithmic and theoretical advances will be 
required in order for such calculations to become a reality. One important ingredient in such 
calculations will be extrapolations of three-nucleon observables as a function of pion mass. 
In particular, a precise understanding of the behavior of three-nucleon observables in the 
critical region will be indispensable in making contact with the physical pion mass. In this 
paper we have investigated the physics involved in such extrapolations and demonstrated 
the complementarity of xEFT and EFT(^) calculations in their development. 

Finally, the possible limit cycle in QCD is also very interesting in its own right. First 
signatures of such a limit cycle have recently been seen in an experiment with cold atoms 
[53j , and it would be interesting to observe a limit cycle in lattice simulations of the three- 
nucleon system with pion masses around 200 MeV [54j . The pionless effective theory can be 
understood as an expansion around this limit cycle and will be instrumental in interpreting 
the lattice results. In this paper, we have demonstrated that the required higher-order 
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calculations are feasible and converge well. 
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